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Abstract. Loop quantum gravity, a non-perturbative and manifestly background free, 
quantum theory of gravity implies that at the kinematical level the spatial geometry is 
discrete in a specific sense. The spirit of background independence also requires a non- 
standard quantum representation of matter. While loop quantization of standard model 
fields has been proposed, detail study of its implications is not yet available. This review 
aims to survey the various efforts and results. 
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1 Introduction 

For most physical situations, we have the comfort of a background space-time geometry - be 
it the space and time of the non-relativistic regime or the Minkowski space-time of special 
relativistic regimes or even more exotic black hole space-times. Even to have a well defined, 
causally well behaved and deterministic classical theory of fields, the background space-times are 
limited to the class of globally hyperbolic geometries. Quantum field theories can be defined in 
all such space-times [31]. However, the familiar Fock representation or the particle interpretation 
can be defined uniquely only for a sub-class of space-times namely those that are stationary or at 
least asymptotically (in future and past) stationary. This is because in relativistic theories, the 
Fock representation depends crucially on having a notion of positive (negative) energy solutions, 
i.e. a notion of preferred time [63]. 

However, in regimes of early universe or evaporating black holes, we loose this comfort of 
having a fixed background space-time. In these cases, not only do we have to face 'quantum 
gravity' but also face the issue of quantization of other matter quantum fields in the absence of 
a background space-time. 

One strategy to construct quantum theories is to use the Gelfand-Naimark-Segal theory [23]. 
Given a commutative C* algebra, for every positive linear functional on it, there is a cyclic, uni- 
tary representation. The strategy is to look for a suitable sub-algebra of the Poisson bracket 
algebra of classical observables, identify its commutative sub-algebra and construct various rep- 
resentation. Choose among these, those on which the full sub-algebra is suitably represented. 
Under favourable conditions, a unique representation is picked out. A well known example of 
this is the unique, weakly continuous representation of the Weyl-Heisenberg C* algebra (non- 
commutative) [4, 19]. The other example is the holonomy-flux representation of Loop Quantum 
Gravity (LQG) [24, 45]. 



*This paper is a contribution to the Special Issue "Loop Quantum Gravity and Cosmology" . The full collection 
is available at http://www.emis.de/journals/SIGMA/LQGC.html 



2 



G. Date and G.M. Hossain 



A natural commutative C* algebra is constructed from suitable configuration space variables. 
For a linear configuration space such as , the coordinates themselves form suitable variables. 
For a more general non- linear configuration space, coordinates have only local meaning and are 
not suitable variables. Instead, one constructs a set of functions on the configuration space 
which are invariant under change of local coordinates. These will be more in number than the 
dimension of the configuration space (manifold) and may possibly have some relations among 
them. For example, for the configuration space S^, we can choose 3 functions on it - namely 
X'^{9,(p), i = 1,2,3 satisfying = 1 (say). The chosen class of functions is required 

to constitute a separating set of functions so that these functions serve to distinguish different 
points of the configuration space. A set of functions {/»} on a space Q forms a separating set 
for Q if, for every pair of distinct points, there are at-least two functions in the set which have 
distinct values at these points. 

For field theories, there is a further feature. The configuration space is heuristically described 
in terms of (tensor etc) fields on a manifold. The 'delta function' in the canonical brackets, means 
that the canonical coordinates are to be understood as suitably smeared fields. For example, for 
scalar field cp, Ff{4>) := d^xf{x)(p{x), for all scalar densities of weight 1, /'s, defines functions 
on the configuration space. For a gauge field Tjj4^dx", the holonomies he{A) := Pexp(/g A), is 
another example of such functions. 

It turns out that suitable representation of the Poisson bracket algebra can also be constructed 
using a special class of functions on the configuration space, the so called cylindrical functions, 
in conjunction with projective techniques [6]. When applied to formulation of gravity in terms 
of the real SU(2) connection, this constructs the kinematical Hilbert space of LQG^. On this 
Hilbert space, geometrical operators have discrete spectra and their eigenstates are associated 
with appropriate class of graphs. This is paraphrased by saying that quantum geometry lives on 
graphs. The quantum excitations of matter are thus also expected to be associated with vertices 
and edges of graphs. 

Construction of such a representation constitutes only a first step in the construction of 
a quantum theory of gravity. Due to the diffeomorphism invariance, the system has gauge 
invariances and the (Dirac) quantization procedure is split into three stages: construction of 
the kinematical Hilbert space followed by (distributional) solutions of the constraints followed 
by specification of physical Hilbert space. 

In this review we will discuss the construction of kinematical Hilbert space. Our discussion 
is confined to four space-time dimensions wherein the matter fields can be scalars, spinors 
and 1-form gauge fields. We will not discuss the gravitational sector (discussed in [6] for 
instance) except to use the real SU(2) variables - connection and densitized triad - when 
needed in the matter sector. Our discussion is within a canonical framework leaving out 
discussion of matter in a spin foam framework. Each type of matter is discussed in a sep- 
arate section. Section 2 discusses loop quantization of scalars and some of its implications. 
Both cases of scalar fields ranging over a compact set (Higgs scalars in the adjoint repre- 
sentation of a gauge group [58]) and those ranging over the full R [7] are discussed. Sec- 
tion 3 describes fermions, their Hamiltonian formulation and loop quantization while Maxwell 
and Yang-Mills fields are discussed in Section 4. Propagation of matter waves on quan- 
tum geometry is briefly discussed in Section 5. Finally, in Section 6 we make a few re- 
marks. 



'^The cylindrical functions in this case are finite linear combinations of spin network functions. To define 
these, take an oriented, piece-wise analytic, closed graph 7, embedded in the 'spatial manifold', E3. Associate to 
each of its E edges an SU{2) group element in a representation j. These are the holonomies of a (generalized) 
connection. Associate to each of its V vertices, an intertwiner c and finally take the trace (contract all indices). 
These functions of connections are naturally labelled by (7,^1, . . . , js, Ci, . . . , cy). Similar functions are defined 
for matter fields as well. 
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2 Scalars 

In loop quantum gravity approach, the scalar fields, like other matter fields, remain compara- 
tively under-studied. Often the scalar field has been used to play a secondary role for describing 
dynamics of gravitational degrees of freedom. In particular, it has been used as a 'clock' variable 
in loop quantum cosmology. Nevertheless, there have been some notable efforts for studying 
quantum scalar fields themselves. In this section, we survey these different efforts involving 
quantum scalar fields. 

2.1 Scalar field as a 'clock' 

We begin the review for scalar fields by considering the situations where scalar fields have been 
used as 'clock' variables. In these studies, mostly for cosmological scenario, scalar fields have 
been quantized using standard quantization. Of course, one should use loop quantization also 
for scalar fields, as done for gravitational degrees of freedom. 

In a general background the equation of motion of a massless, free scalar field cp is 

d^{V^g^''d,^)=0, (2.1) 

where g is the determinant of the metric ^^jy. For cosmological backgrounds, i.e. in a spatially 
homogeneous space-time, the equation of motion (2.1) reduces to 

^t{^/^N~^^t(|)) =0. (2.2) 

Now if one makes the choice of lapse function N to be \/q{t) where q is the determinant of the 
spatial metric then the equation (2.2) implies scalar field (j) (x t. This suggests that one may 
choose a massless scalar field as an 'internal clock' to describe the dynamics of the remaining 
degrees of freedom. Use of 'internal time' to describe the dynamics, instead of the coordinate 
time t, has the advantage of avoiding the so-called 'problem of time' in quantum gravity. 

By choosing such a lapse function, the total Hamiltonian density for a cosmological system 
with a massless free scalar field can be written as (A^ = ^/q chosen) 

H = Hg+ ]^TT% 

where Hg represents gravitational Hamiltonian. Now imposition of Hamiltonian constraint on 
the Hilbert space, i.e. = implies 

= 2Hg^, 

where we have used tt^ = —id^. Thus, the positive and negative frequency solutions can be 
viewed as satisfying a set of first order evolution equations of the form 



By identifying the field 4> to play the role of time and IH that of a proper Hamiltonian, one may 
note that the positive or negative frequency solutions satisfy appropriate Schrodinger equation 
with respect to the internal time (j)- This approach of considering scalar field as an internal 
time has been used to describe the quantum dynamics of isotropic cosmological models such as 
Friedmann-Robertson- Walker space-time [8] as well as anisotropic cosmology such as Bianchi I 
space-time [11]. It has been shown for isotropic models that universe undergoes a bounce [9] 
when quantum dynamics is described with respect to an internal time. 



4 



G. Date and G.M. Hossain 



2.2 Compact scalars 

The Higgs scalars plays a very important role in standard model of particle physics. It provides 
a mechanism for generation of masses for different elementary particles. Naturally, in any 
approach to quantum gravity, one would expect to have an appropriate representation of the 
Higgs scalars in the given quantum theory. 

In the standard model of particle physics, while constructing the quantum framework one 
specifically uses the fact that the background space-time is a Minkowski space-time. A natural 
kinematical measure available in constructing a suitable Hilbert space in a such space-time is 
Gaussian measure. Use of this measure leads to the standard Fock Hilbert space. However, one 
cannot use Gaussian measure in a background independent construction of the Hilbert space such 
as those used in loop quantum gravity. This is because Gaussian measure is not diffeomorphism 
invariant. 

In order to avoid this problem, Thiemann proposes [58, 59, 60, 61] to use a set of bounded 
variables for quantizing Higgs scalars. In particular, bounded variables of the form U{x) = 
^(I>i{x)t^ are used as the configuration space variables rather than the scalars (/>/ themselves. 
Here, scalars <f>i have been made dimensionless with the help of a constant scale of VSttG where 
G is the Newton's constant of gravitation. These variables are the analogues of the holonomies 
of connection variables. Unlike the gauge holonomies which are labelled by curves, these are 
labelled by points (no smearing) and are referred to as 'point-holonomies'. These are valued 
in Q since the 4>j{x)t^ are valued in its Lie algebra. Note that U{x) is a matrix representing 
a group element in a certain representation determined by the generators . Matrix elements 
of U{x) therefore determine functions of the Higgs scalars. If one considers a single real scalar 
field, valued in a compact interval [0,27r] (say), then the basic variables can be chosen to be 
U{x) = e*<^(^') which are valued in U{1). For constructing the kinematical Hilbert space, one 
considers products of finitely many such matrix elements labelled by a finite set of points of 
the spatial manifold together with unitary representations of G associated with them. Finite 
linear combinations of these functions generate the (quantum) configuration space U and the 
kinematical Hilbert space is then obtained as the space of square integrable functions L2{V(, d^n) 
where dfiH is naturally available Haar measure on the group Q. On this, the point holonomies, 
but not the scalar field, are well defined (multiplicative) operators. The momentum is a density 
weight 1 object and its integral over three-dimensional regions are analogues of the flux variables 
of the gravity sector. Their action on the polynomials needs careful definition for which we refer 
to [58]. All other operators of interest are to be constructed using these basic operators. 

We may recall that the classical dynamics of a real scalar field is governed by the Hamiltonian 



Ha = / d^xN 



(2.3) 



Formally, one can construct the quantum Hamiltonian operator corresponding to the classical 
Hamiltonian (2.3). Regularization of this quantum Hamiltonian operator is performed by using 
the triangulation of manifold which is analogous to the construction of gravity sector opera- 
tors. While formal construction of the scalar Hamiltonian has been performed by Thiemann, 
understanding it's physical consequences remains an open issue. 



2.3 Non-compact scalars 

In Thiemann's construction [58] , the configuration space variables are essentially periodic func- 
tions of the scalar field. This is sufficient to describe the scalar fields which range over a compact 
interval. However, such treatment of scalars is inadequate to describe the Klein-Gordon type 
real- valued scalar field which takes values in rather than in U{1). To circumvent this prob- 
lem, Ashtekar, Lewandowski and Sahlmann propose the use of almost periodic functions as the 
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configuration space variables [7]. The formal construction of the kinematical Hilbert space for 
a single non-compact scalar is carried out as follows: First one considers finite sets of points 
in referred to as the vertex sets. For a given vertex set V = {xi, . . . , the corresponding 
vector space is generated by linear combinations of the functions of the form 

where Aj's are real numbers. After considering all possible sets of vertices, the union of 
corresponding vector spaces and together with its completion, one construct the quantum 
configuration space required for loop representation. If one denotes this space as A then the 
corresponding Hilbert space is given by L2{A,dfi). The measure is defined as follows 

d,«,,= ^- «^^=»"^- (2.4) 
A ' I 0, otherwise. 

The elementary configuration operators are analogous of the holonomy operator for connection 
variables and they act by multiplications as 

Use of the measure (2.4), makes the elementary holonomy operator h{x,X) = e'^^'^^^^ not to 
be weakly continuous in the parameter A. This implies that one cannot construct an operator 
corresponding to the scalar field </>(x) itself. The conjugate variable to the holonomy which is 
promoted as an operator in the quantum theory, is taken to be the smeared field momentum 



Pif) = J dx-'TT^{x)f{x), 

where vr^ is conjugate momentum field and f{x) is a test function. Their Poisson bracket is 
given by 

{h{x,X),P{f)}=i\f{x)h[x,\). (2.5) 

Subsequently, one looks for a representation of the Poisson bracket (2.5) in the quantum theo- 
ry as a commutator of the elementary operators. The non-trivial commutator bracket of the 
elementary operators is 

[h{x,\),P{f)]=iXf{x)h{x,X), 

as holonomy operators and smeared momentum operators commute with themselves. After 
having the kinematical framework ready, one may construct the scalar Hamiltonian operator to 
understand its quantum dynamics. 

2.4 Parametrized scalar fields 

Within loop approach, scalar fields have been studied in the context of parametrized field theory 
(PFT). In the PFT approach, scalar field theory on a flat space-time is recast in the form as 
if the fields reside in a generally covariant background. In the references [43, 44], Laddha and 
Varadarajan have studied parametrized field theory for a free scalar field cj) in 2-dimensional 
flat space-time. The action for such a scalar field with respect to an inertial frame is 

S[4>] = [ d^Xi^^^'dAcpdBcl), (2.6) 
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where metric ijab = diag(— 1,1) and A = 0,1. Now one may choose to use an arbitrary 
coordinate system involving x"" with a = 0, 1 such that can be viewed as 'parametrized' 
by x"' as = X'^{x°'). Then the action (2.6) can be rewritten as 



S[(t>] = j Sx^g'^^'daCt^da^ 



where gab = VabdaX^dbX^ and g is the absolute value of the determinant of the metric gab- 
In the parametrized field theory, the action £'[(/>] is treated as not only a functional of (p, but 
also a functional of X^^s. In other words the action for parametrized field theory is taken to be 

S[(t>,X^] = [ (fx,/^)g''\X)da<Pda<p. (2.7) 



The variation of the action (2.7) with respect to the field 4> leads to da{^g°'''4>da4>) = 0. This 
equation is basically the same as the flat space equation of motion for cj), i.e. r]^^dAdB4> = but 
written in coordinates. The variation with respect to X^ leads to the equations which are 
automatically satisfied as long as r]^^ dAds'P = 0. In other words, choice of the new coordinates 
remains undetermined as expected. 

For canonical formulation of this theory, it is convenient to choose the light-cone coordinates 
defined a.s X^ = X^ ± X^ . One may denote the corresponding conjugate momenta by II-i-. In 
terms of these coordinates the action (2.7) can be expressed as 



s = j dt j dx[i:^4) + n+x+ + n^x- - n+h+ - n~h^ 



where A^^ are the Lagrange multipliers for the respective constraints H±. As constructed, the 
non-trivial Poisson brackets are 

{(/.(x), 7r<^(x')} = 5{x, x'), {X±(x), n±(x')} = 5{x, x'). (2.8) 

In the standard Fock quantization, one seeks an appropriate representation of the Poisson brac- 
kets (2.8) as commutator brackets of the operators in the quantum theory. 

To proceed to polymer quantization, one notes that in 1-space dimension any tensor can be 
viewed as a density of weight equal to its covariant rank minus its contravariant rank. Conversely, 
a density weight 1 object can be viewed as a 1-form. Hence II-i-, Y± := ^ dx4> are all 1-forms 
which are integrated along intervals. Exponentials of the integrals are the holonomies {not point 
holonomies). The polymer quantization of the system consists of choosing holonomies of II-i- 
as multiplicative operators for the X^, Il± sector and choosing holonomies of both Y± in the 
matter sector as the basic variables (non-commuting). Details should be seen in the references 
cited above. 

Using this quantization it is shown that one can construct a suitable state in the polymer 
Hilbert space which reproduces the Fock space two points function in the large wavelength 
limit. However, understanding the behaviour of two-point function in the strong polymer regime 
remains an open issue. 



2.5 Polymer quantized scalars in Fourier space 

A different study of scalar fields using loop quantum gravity techniques has done by Hossain, 
Husain and Seahra [34]. In their approach, polymer quantization is performed for the Fourier 
modes of the free scalar field in flat space-time. This has the advantage of dealing with a known 
Hamiltonian which is a sum of decoupled harmonic oscillators. One then computes 2-point 
correlation function directly in the energy eigenfunction basis. 
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The phase space variables for the free scalar field are the canonical pair ((/)(x, t), 7r(x, t)) 
satisfying the Poisson bracket 

{</.(t,x),7r(t,y)} = <5(3)(x-y). 
The Hamiltonian is 



(2.9) 



^ 1 ab 



2 +2^^5aM 

where the space-time metric is ds"^ = —dt^ + rjabdx"'dx^ . It is convenient though not essential 
to put the system in a box by restricting to a finite region of flat 3-space so that the volume 
V = j d^Xy^ is finite. One can then perform Fourier expansion of the field into the 3-momentum 
space as 



<A(t,x) = ^ J]</,k(^)e*'^■^ Mt) = [d'x e-''-(/.(x,t), 



along with a similar expansion for 7r(x, t). After a suitable redefinition of the Fourier modes, 
one can express scalar field Hamiltonian (2.9) in terms of real-valued (j)]^ and vr^ as 



k 



with the Poisson bracket being {(/'k, '?i"k'} = i^kk'- Clearly, the Hamiltonian is a sum of decoupled 
harmonic oscillators. The polymer quantization of each mode then follows that of a quantum 
oscillator. In particular, the variables used in polymer quantization are (/)k and f/Ak = e*^'^'= 
which satisfy the Poisson bracket {^k, f^Ak} = i^Uxu.- The parameter A has dimensions of 
(length) 1/2. 

As in the case of polymer quantization, vr^ itself cannot be a well-defined operator in the 
quantum theory as the action of C/^k is not weakly continuous with respect to A. To represent 
momentum operator, instead one uses vr^ = (C/A*k — Uj^ k)/^*'^*- emphasized that the 

polymer quantization method introduces a new length scale (which is A^ here) in addition to 
Planck's constant into the quantum theory. One can define a dimensionless parameter as 

2 1, I |k| polymer length scale 

g = XJk\ = -T-r fx — — , 

spatial wavelength 

where is the fundamental length scale associated with the polymer quantization of (p. 

Clearly, g = should recover the results that one get from Fock quantization. The standard 
2-point correlation function is defined as 



(O|<A(x,t)0(x',t')|O) = :^5]e^''-("-"')Z)k(t-t' 



V 

k 

where |0) = Hk |0k) is the vacuum state and the matrix element is 

D^{t - t') = (Ok|e^-^''*0ke"*-^''*e*^''*'4e"^-^''*'|Ok), (2.10) 

where is the Hamiltonian operator. The matrix element (2.10) can be computed using the 
polymer oscillator spectrum -ffkl^k) = E'^\n\f^, and the expansion of the state 0k I Ok) in the 
energy eigenstates as </>k|Ok) = X]n'^"l'^k)- By defining 4-momentum p = (a;,k), one can write 
the momentum space propagator as 

p2 + AS2_|k|2_^g' (2.11) 
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where p = —00 + |k| and choice of the sign of ie corresponds to the Feynman propagator. One 
can recover the expected Fock space result Dp = i/{p'^—ie) by using properties of the Schrodinger 



oscillator, i.e. Cn = 2|k| and AEn = n|k|. The propagator (2.11) with polymer corrections 

in the infrared limit (5 <C 1) can be written as 

— 5|k|^ — le ' 
The pole of the propagator implies an effective dispersion relation 
io^ = |k|2(l - \k\/AQ, 

which violates Lorentz symmetry. On the other hand in the ultraviolet limit {g ^ 1), the 
propagator with leading order corrections is 

D - ^^^^"^ O ( ^ 

^ p2 _|_ 4^2|J,|2 _ \^^6 



The corresponding dispersion relation 
= 4\k\^/M^ 

also violates Lorentz invariance. One may also note that the propagation amplitude at high 
momentum is suppressed by a factor 1/g^. 

Polymer quantization introduces a length scale = M~^. This is unavoidable. Such a scale 
could be provided by an underlying quantum geometry. However at the present level of under- 
standing, it is to be treated as an arbitrary parameter. Current observational limits on Lorentz 
violation may be then used to constrain its value. 



2.6 Other studies of scalar fields 

Ultraviolet behaviour of polymer quantized scalar fields has been studied by Husain and Kreien- 
buehl by constructing Fock-like states in flat space-time [35]. Instead of Fourier space, as dis- 
cussed in section (2.5), here one applies polymer quantization in real space. It is shown that the 
vacuum expectation values of the commutator and anti-commutator of the creation and annihi- 
lation operators become energy dependent, and tends to show some sort of fermionic behavior 
at high energy. Furthermore, the modified dispersion relation that arises leads to violation of 
Lorentz invariance. 

In the context of spherically symmetry, Gambini, Pullin, and Rastgoo have studied polymer 
quantization of scalar field [27] coupled to gravity. Their study is performed in the midi- 
superspace context where even after imposition of the given symmetry, field theoretical nature of 
scalar fields, i.e. with infinitely many degrees of freedom, is retained. After choosing a quantum 
geometry state corresponding to locally flat geometry, quantization of the spherically symmetric 
scalar field is carried out using a discretization followed by a 'polymerization' done in two 
different ways (the field and the momentum). Their notable conclusion is that the propagator 
obtained is not Lorentz invariant. 

Polymer quantized free massless scalar field in a homogeneous and isotropic cosmological 
space-time has been studied by Hossain, Husain and Seahra [33]. Use of semi-classical Friedman 
equation yields a non-singular and non-bouncing universe, without quantum gravity. The system 
exhibits an early de Sitter-like inflationary phase with sufficient expansion to resolve the horizon 
and entropy problems, and a built in mechanism for a graceful exit from inflation. 

The propagation of polymer quantized scalar field in flat space-times has been studied in [32] . 
Polymer quantization is performed in real space there unlike in Fourier space as discussed in 
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section (2.5). Using semi-classical states, it has been shown that effective wave equation is 
both nonlinear and Lorentz invariance violating. Furthermore, it is demonstrated that polymer 
effects tend to accumulate with time for plane-symmetric waveforms. Possibility of measuring 
deviations from the Klein-Gordon equation in particle accelerators or astrophysical observations 
is also discussed. 

3 Fermions 

Fermions coupled to gravity has been discussed in the literature [12, 29, 39, 49, 50, 51, 64], at 
both classical and quantum level. In the formulation of general relativity in terms of real SU(2) 
connection, Thiemann discussed loop quantization of standard model fields [58, 59, 60, 61]. 
The fermions were treated in the second order form, i.e. fermions couple to gravity through the 
spin connection (torsion free Lorentz connection). Perez and Rovelli returned to fermions in 
presence of the Hoist term and found that the Barbero-Immirzi parameter, 7, inverse of the 
coefficient of the Hoist term, becomes classically observable [53] . Mercuri [47, 48] discovered that 
with a further addition of suitable non-minimal fermionic couplings, 7 can be made classically 
unobservable. He also noted that the added terms (Hoist plus non- minimal) can be expressed as 
the Nieh-Yan topological term after using the connection equation of motion. This strategy of 
adding non- minimal couplings to keep 7 classically unobservable was followed for = 1 , 2 and 4 
supergravities also [41]. Canonical analysis and loop quantization of fermions with non-minimal 
couplings was discussed by Bojowald and Das [14, 15]. 

It was subsequently realised that 7 will automatically be classically unobservable provided it 
is the (inverse) of the coefficient of the Nieh-Yan term (a total divergence) in the Lagrangian 
density. Thus, instead of the Hoist terms alone, if the Nieh-Yan term (Hoist + (torsion)^ piece) 
is used in conjunction with the Hilbert-Palatini, then for arbitrary matter and their couplings, 
the 7 will drop out of the classical equations of motion^. But now, since the action is modified, 
it was not obvious that the real SU(2) formulation will result from the new action. It turned 
out that it is possible to systematically derive the real SU(2) Hamiltonian formulation from 
such an action [21, 57]. Since it is in presence of fermions that non-trivial torsion results from 
the equation of motion of the Lorentz connection, fermions were also included in the canonical 
analysis and real SU(2) formulation was seen to emerge. The canonical analysis leading to real 
SU(2) formulation has since been extended to include the other two topological terms namely 
the Pontryagin and the Euler classes [42, 54]. 

It is straight forward to derive the real SU(2) formulation from the Hilbert-Palatini action 
with 7~^ times the Nieh-Yan term [20, 21] . Its main points may be summarised as follows. In the 
time gauge parametrization^ , to begin with there are 13 components of the co-tetrad (A^, A^", V^) 
and 24 components of the Lorentz connection (w/'^, uj^'^)- Of these, 10 variables - N, N""^ ajj."^ - 
occur as Lagrange multipliers and 18 variables - Ef ~ Vf" (densitized triad) and 9 combinations 

of u^"^ - explicitly appear in the form pq. The coefficients of the Lagrange multipliers form 
the primary constraints Tia, Ti, and Q"^^ . The remaining 9 combinations of the Lorentz 
connections, have no velocities and lead to additional 9 primary constraints, say vTj ~ ~ vr*-^ 
(symmetric in i, j). Of these, the 3 primary constraints, vTj ~ form a second class system with 
the primary constraints Q^^ and are eliminated easily using Dirac brackets. Preservation of the 
remaining 6 primary constraints vr*-' w lead to further 6 secondary constraints Sij ~ 0, with 
which they form a second class system. These are again eliminated using Dirac brackets. The 
left over system has the 18 phase space variables, (AJj,£'?), and 7 first class constraints and is 

necessary condition for a topological origin of 7 is thus satisfied. 

J = 0, 1, 2, 3 are the Lorentz indices, a — 1,2,3 denote the spacial indices. The co-tetrad components are 
et := A'^, e] := N'^Va, ej, := V^, qab '■= V^V^Sij, q :— det{qab) 7^ 0. The Lorentz connection components are 
denoted as 1^°' =: Kl, oj'J =: E'^V^ [20]. 
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the real SU(2) connection formulation. This counting and the steps in the constraint analysis 
remain the same when matter is included. 

This derivation also introduces factors of sgn(e) := sign(det(e^)) (= A^sign(det(14*)) in the 
time gauge parametrization) in appropriate places. For instance, we are naturally lead to the 
definitions: := sgn(e)y^V^" and A^^ := 7~^sgn(e) -fC* — TJ^ Under an improper orthogonal 
transformation A*^ G 0(3) acting on the index i, the triad changes its handedness and the 
sgn(e) factor changes sign leaving the handedness of Ef unchanged. This is as it should be 
since the index i on E represents adjoint representation of SO (3) while on V it represents the 
defining representation. For S0(3), both are equivalent but not for 0(3). Under an inversion, 
A^j = quantities in the defining representation change sign while those in the adjoint 

don't. The same reasoning applies to the definition of the connection. Now the connection is 
also even under inversion. The sgn(e) factor also change the behaviour of Ef and under 
the action of orientation reversing dif feomorphisms. These factors are relevant for discussion of 
'parity' properties of the canonical formulation (discussed later). This is independent of matter 
couplings. 

When a Dirac fermion is included, the solution of the secondary, second class constraint 
Sij ~ changes since fermions couple to Lorentz connection apart from coupling to the triad. 
This solution leads to non-trivial Dirac brackets between the SU(2) connection and the fermions. 
One can however make natural shifts in the definition of the connection to recover the canonical 
brackets. This also simplifies the constraints. Four Fermi interaction terms however survive in 
the Hamiltonian and are signatures of first order formulation. In the second order formulation 
where fermions couple to the torsion free connection, there are no terms quartic in the fermions. 

Fermions are also tied with possible parity violations [14, 15, 25]. There are two distinct 
notions of 'parity': one related to orientation of the space-time manifold (parity) and one related 
to the improper Lorentz transformation (Lorentz parity). Depending upon the definitions of the 
basic canonical variables (with or without the sgn(e) factors in this work (Section 3.3)), the 
canonical framework and the action are (non-) invariant under one of the notions of parity. 
These possibilities are distinguished and discussed in [20]. Bojowald and Das discuss the non- 
invariance under Lorentz parity [14, 15] in the context of non-minimally coupled fermion. 

We now turn to the classical, canonical form of a Dirac fermion minimally coupled to gravity 
in the first order formulation and discuss the loop quantization of fermions in the following 
sub-section. 

3.1 The Hamiltonian formulation 

The starting point is a choice of (tensor /spinor) fields and a corresponding generally covariant, 
local action on 4-dimensional space-time M ~ R x S3. The next step is to carry out a 3+1 
decomposition to identify the Lagrangian which is a function of (tensor) fields on S3 together 
with their velocities with respect to the chosen time coordinate. The fields whose velocities 
appear in the Lagrangian are potentially the configuration space variables while those with- 
out velocities appearing in the Lagrangian are Lagrange multipliers whose coefficients will be 
primary constraints. This Lagrangian leads to the kinematical phase space. Now a constraint 
analysis a la Dirac is performed. If there are second class constraints, one may hope to simplify 
the analysis by solving the second class constraints. However, now one must use the Dirac 
brackets. These may not have the canonical form for the remaining variables (i.e. may not be 
Darboux coordinates) and a new choice of variables may be necessary. This is particularly rel- 
evant for Lagrangians which are linear in the velocities eg the Hilbert-Palatini-Nieh-Yan and 
the Dirac Lagrangians which typically do have second class constraints. The classical Hamilto- 
nian formulation is completed when the action is expressed in the Hamiltonian form together 
with first class constraints. We also have fields coordinatizing the kinematical phase space 



Matter in Loop Quantum Gravity 



11 



(after the second class constraints are eliminated) with the configuration space coordinates 
identified. 

We begin with the Lagrangian 4-forms built from the basic fields the co-tetrad e^^dx^ and 
the Lorentz connection ujj/dx'^ {k := SirG) 



2k 



£NY(e,a;) = [T^ {e,uj) ATi{e,oj) - Rij{oj) A A e'^] 



-^Dirac — 2 I ^ 



■ ~Xe';j'D^{oj,A,...)X-D^{co,A,...)Xe'',yX 
1 

1 T J i ./.,+ . .,+1 ^ 

7 • 



D^{uj)X := d^X + ^uj^^aijX + ie'A^^X + • • • + A, 
D;M^ := ja^At + ^^-^Ata]^ - ie'A^X^ + . . . + Atj 7° 



The . . . refer to possible couplings of the Dirac fermion to other gauge fields, e.g. the Maxwell 
field. These are suppressed in the following'^. 

Here R := du + uj A u and T := de + u A e are the usual curvature and torsion 2-forms. The 
factor of sgn(e) is present because only then the Hilbert-Palatini Lagrangian matches with the 
y^\g\R{g). This arises from noting that determinant of the co-tetrad is given by, e = sgn(e) y^H- 

A 3+1 decomposition is carried out as usual by choosing a foliation defined by a time function, 
T : M — )• [R and a vector field f^d^, transversal to its leaves. The vector field is normalised by 
t ■ dT = 1 so that the parameters of its integral curves, serve as the time coordinate. Given 
such a decomposition, we choose a parametrization of the co-tetrad which leads to the usual 
ADM parametrization in terms of the metric, gfj_,y := ejj^elrjjj. Both the co-tetrad and the 
corresponding tetrad in this parametrization are displayed below 



ei = Nn' + N^V^, ei = Vl n'ni = -\, n^Vj = 0, (3.1) 
with VfVi = Sj + nin\ V^V^ = 8^. (3.2) 



We will now restrict to configurations such that rii = 0, no = —1. This also implies that 
V^'' = = and that are invertible with as the inverse^. We also define the 3-metric 
(lab '■= ^d^b^ij (which is positive definite in classical theory) and denote q := det{qab)- 

For future convenience we introduce \I' := q^^^X, := g^^^A^. This absorbs away the 
factors in the Lagrangian as well as in the constraints. Note that the terms involving 
the derivatives of ^/q cancel out. The A fermionic variables being of density weight zero, the ^ 
fermionic variables are of density weight 1/2. From now on we will use the half density variables. 

Substituting the 3+1 parametrization of the tetrad and using the time-gauge, the Lagrangian 
can be written as 

CDir.c = li^^dt^ - dti^^m - u;mG'P - \u^tij£''''G^ - N'^'n^, - NTip, 

''Our conventions are: Greek letters denote space-time indices, lower case Latin letters denote space indices, 
upper case Latin ones denote Lorentz indices. Our signature is ( — h H — h). The (metric independent) Levi-Civita 
symbols are S'^^^' = +1 = £:oi23. The determinant is defined by e f^J^^ = -^'"'"''e^ei^ef e^. For the Dirac 
matrices: 2ri"l = j'Y + 1'' ^' , ■= 3 [7^,7'^], 75 := i7°7S^7^, A := A'l"7''. 

^This would correspond to the choice of the so-called time gauge if we started without restricting the 
configurations a priori. 
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where 



niF 

rL„i 



For the fermions, the action being hnear in velocities, we have primary constraints, vr^ ~ A^, 
7r;^t ~ A, which are second class. Also these variables fail to be Darboux coordinates ~ do 
not have vanishing Poisson brackets^ with the gravitational variables due to the factor. The 
shift to ^, variables makes the matter and gravitational variables Poisson-commute. Defining 
Dirac brackets relative to these primary, second class constraints allows us to use ^'^ as basic 
variables with Dirac brackets given by 

From the details given in [20], we have the final expressions: 

P,- := {K^r'Et = iKj)-hgn{e)V:^^, = 7sgn(e)i^^ - Ti{V), 

fijk ^ 



-yk{dbVaj - daVb, + VfV^dbV^i}, 



. 2 - ^J>1A{£^)fUA) - (1 + 7^) {KiK^, - KiKt) } + jda (sgn(e)y«g;,,) 



v. := K7 



+ 



2 

1 P^P^ 
2 

i K7sgn(e)Pj" 



^/9 



3_ (^t^5^) 
16" ^ 



{^YVa^ -Va^f^}{A) + 
2 ($t^5^0i^)(^t^^^0.^y 



K7 



K'P'' 



^'"^75^) 



(3.3) 



In the above, = OaP,^ + S.-AiP^ := - g>. 



Remark 1. Dimensionally, k ~ L^, {A,K,d) ~ L~^, E ^ , P ^ L''^, {^,^) ~ L'^^'^, 

The density weights are: {P,g,'Ha,'H) = +1, (^,*) = 1/2 and {A,K,V,r) = 0. 

Under'' Lorentz parity: {V, K, sgn(e)) are odd, (F, A, P, ^) and the constraints are all even. 

Under parity combined with 7 — )• — j: {V, K, F, A, P, 7 sgn(e)) and the constraints are all even. 

The A and K above correspond to the vacuum case for which the Thiemann identities hold. 
The inverse square root of q and i^* appearing above are manipulated exactly as in the vacuum 
case. 

Explicitly, the identities are 



sgn(e)f^-K 
sgn(e) 



gijk_ ^j^k 



K7- 



Vd^tpf)' 
= <i K{x), [ d'y^ 



q := {det{V:)r = det{Ef 



^Strictly Generalized Poisson brackets [30], due to the Grassmann nature of the fermions. 
'^Discussed in Subsection 3.3. 



Matter in Loop Quantum Gravity 



13 



E E'^ 1 ( /" 1 

S'^^ J ^ = —S^^" { 4 (a;), / d^yJ^ \ , 

He{1) := / ^^'X- ^ -.= 1 d'ysgn{e)KPt 

K = {kj^)-' |i?E(l),y" d'y^^ , sgn{e)Klix) = {Al{x),K}. 

These identities suffice to derive a quantization the Hamihonian constraint from that of the 
'Euchdcan Hamihonian constraint' (the first term in the Hamihonian constraint) and of the 
volume operator. This completes the classical canonical formulations as it follows from the 
action. 



3.2 Constraint algebra 

It is easy to see that the gauge constraint generates correct gauge transformation of the basic 
fields. Specifically, with ^(A) := /j,^ d^xK^Qi, 



{A\{x),gm = -T^aA' = -daA'-S'jkAiA", 
{P^{x)MA)] = +£^fA=P^, 
{*"(x),g(A)} = -iA*(75CTo^^r, 
{^i{x),g{A)} = +iA\^h^<y^i)a- 

If we compute the infinitesimal action of the Ha constraint on the basic variables, we see that 
it equals the Lie derivatives of the basic variables only up to an SU(2) gauge transformation. 
We are however free to modify the constraints by adding suitable combinations of themselves. 
So we define the diffeomorphism constraint as 



C{N) := [ . 



d^xN'^Ca with 



Ca := Ha - AlGi = P^'daAl - dbiA^P^) + '-{^^da^ - Oa^^ ■ *), 

which leads to the infinitesimal transformations, 

{Aj,(x),C(iV)} = C^K = 9a(iV'^l) + N\dkK - da^,), 
{P^{x),C{N)} = Cf^P^ = N%P^- - P^d.N'^ + 1 ■ {dbN>^)Pt, 

{*"(x),C(7V)} = Cj^^'^ = AT^Sft*" + ^ • {dbN'')^'^, 
{^iix),CiN)} = C^^i = N%^i + ^ • id,N')^l 

This implies that {var, / N"'Ca} = Cm^vaj) for all variables. The Gauge constraint already 
generates the correct gauge transformation of the basic variables. By inspection, it follows 
that the gauge constraints (weakly) commute with the diffeomorphism and the Hamiltonian 
constraint, the gauge constraint and the diffeomorphism constraints form sub-algebras and the 
diffeomorphism constraint transforms the Hamiltonian constraint by the Lie derivative. The 
non-trivial bracket is the bracket of two Hamiltonian constraints. 

Before turning to loop quantization, we briefly draw attention to the invariance (or lack of 
it) under two distinct 'parity' operations. 
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3.3 Parity and Lorentz parity 

Recall that we begin with the (co-)tetrad field e^, the Lorentz connection ioj/ and the fermion 
fields A, A (or ^) defined over a manifold M ~ R x S3 which is assumed to be orientable. 
With the topology specified, M can be taken to be time- orientable with respect to all the metric 
tensors constructed by the parametrization (3.1), (3.2). Obviously, II3 is orientable as well. 

There are two distinct sets of 'parity' transformations: orientation reversing diffeomorphism 
of M and a 0(1,3) transformation with determinant = —1. Note that these distinctions do 
not exist in flat space-time with Minkowski metric where the allowed diffeomorphisms are the 
isometrics of the Minkowski metric which are the same as 0(1,3) transformations. We will keep 
the time orientation fixed. Orientation reversing diffeomorphisms of M will then be reversing 
the orientation of S3. We will refer to these as parity transformations. The improper Lorentz 
transformations Aj, will also be taken so that det A = — 1 and A'^q = 1 and will be referred to 
as Lorentz parity transformations. 

After going to the canonical framework in the 'time gauge', we have the fields A^, P^, 
^ defined on S3. A parity transformation is now an orientation reversing diffeomorphism 
of S3 while the improper 0(3) transformation, inversion, will be taken as the Lorentz parity 
transformation. 

In the Lagrangian framework, the Hilbert-Palatini action is invariant under both sets of 
transformations while the Nieh-Yan (and the Pontryagin) actions change signs under parity but 
are invariant under Lorentz parity. Hence the combined action is invariant under Lorentz parity 
and non-invariant under parity. The sgn(e) factor in the Hilbert-Palatini action is crucial for 
this. 

The variables of the canonical framework are defined in terms of those of the Lagrangian 
framework. These definitions of the SU(2) connection in terms of K and F and the conjugate 
momentum in terms of the triad are consistent with the SO (3) gauge transformation extended 
to include the Lorentz parity. Thus the triad which transforms by the defining representation 
changes sign under Lorentz parity. The 'densitised triad' (or the conjugate momentum) trans- 
forms by the adjoint representation and should be invariant under Lorentz parity. The sgn(e) 
factor in their definitions precisely takes care of this. The same can be seen in the definition of 
the connection. It is easy to see that the symplectic structure and the constraints (vacuum) are 
all invariant under Lorentz parity. 

When fermions are included, these are scalars under orientation preserving diffeomorphisms^, 
and transform as ^ — )• 7"^ under Lorentz parity. All the constraints including fermions are in- 
variant under Lorentz parity. This is true in both the Lagrangian and the canonical frameworks. 

With regards to parity the situation is different. The action is not invariant under parity, due 
to the Nieh-Yan term. In the canonical framework, the connection is not simply even/odd under 
parity since the K term changes sign while the T does not. The 'densitised triad' also acquires 
an extra minus sign under parity (behaves as a 'pseudo- vector of weight 1'). The symplectic 
structure thus is not invariant. The constraints also are not invariant under parity and this is 
consistent with the non-invariance of the action. 

The action is invariant under parity combined with 7 — )• —7. Our definitions have the 
appropriate factors of 7 to restore the simple (even) behaviour of the basic canonical variable 
resulting also in the invariance of the Poisson brackets and constraints. 

In short, Lorentz parity is an invariance of the action as well as the canonical framework 
while parity is not. However, parity combined with 7 — t- —7 is an invariance of both action and 
the canonical framework. 

^Definition of spinors depends on the orientation. Thus the action of parity on spinors needs a much more 
careful statement. We will assume for simplicity that the under orientation reversing diffeomorphism, the fermionic 
action remains invariant. 
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Remark 2. We could consider a canonical formulation ab initio, say by making a canonical 
transformation from the ADM variables, without any reference to an action. We could then 
define of the basic canonical variables A^, i??, without the factors of sgn(e). This will restore 
the 'densitized triad' to its usual density weight 1 vector density status and the connection to 
its 1-form status. The canonical framework is then invariant under parity (without changing 
sign of 7). Under Lorentz parity, however, the connection does not have simple behaviour and 
the conjugate momentum is odd. The reason of course is that the internal index 'i' on the triad 
and on the extrinsic curvature refers to the defining representation while that on the connection 
and its conjugate refer to the adjoint representation and these are distinguished by the improper 
orthogonal transformations. This results in non-invariance of the canonical framework under 
Lorentz parity. If the sign of 7 is changed along with the Lorentz parity transformation, then 
the basic variables are even, the symplectic structure is invariant and so are the constraints and 
the action. Thus, the definitions without the sgn(e) factors, interchanges the role of Lorentz 
parity and parity appropriately combined with 7 — t- —7. 

Which of these notions is 'appropriate'? If we were to consider formulation in terms of the 
metric tensor, then the notion of Lorentz parity is not even definable as there is no internal 
Lorentz transformation. On the other hand, with fermionic matter we have to consider ori- 
ented manifolds and the tetrad formulation, introducing the possibility of Lorentz parity. If 
the orientation of the manifold is regarded as a fixed background structure, the parity transfor- 
mations are excluded by definition and Lorentz parity alone remains. Which of these notions 
is relevant for experimental observation, is unclear and so is the issue of 'gravitational parity 
violation/invariance' in LQG. 

In the next subsection, we briefly summarise the loop quantization of fermions [14, 15, 58, 
59, 60, 61]. 

3.4 Loop quantization of fermions 

The loop quantization of the gravitational fields as well as of the fermions has been already given 
by Thiemann and those arguments remain valid. We have already gone from A fermions to the 
half density ^ fermions. In our derivation from the action, this is simply seen as the choice 
of Darboux coordinates. We also do not work with the 'conjugate' variables vrij,, irq,. Instead 
we solve the primary second class constraints and use Dirac brackets relative to these. Thus 
we take, for a single Dirac fermion, "if^, to be complex, Grassmann valued, half densities: 
'■= (^'^)*, A = 1, ... ,4. This can easily be generalized to many Dirac fermions or Weyl 



As mentioned in the introduction, the first step is to construct functions on the fermionic 
phase space. Here we discuss Thiemann's proposal presented in a slightly different form. It is 
based on a particular, naturally available, class of smearing densities. 

Consider a decomposition of S3 into cells Bn, closed subsets of S3. Let := iJ,{Bn) be the 
Lebesgue measure ('coordinate volume') of the cell Bn and let Vn be the centre of Bn (a marked 
point in the interior of Bn)- Define the indicator or (characteristic) function: 



Then the quantity, Xn{x)/e^ has density weight w. This is because, the Lebesgue measure 
of a cell, transforms as a scalar density of weight —1 under a coordinate transformation and 
the indicator function is of course invariant. Note that eventually we are interested in the 

^The overbar on a with an index A, denotes the complex conjugate while one without an index denotes the 
Dirac conjugate: ^ = vpt^". 
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limit of infinite refinement of the cell decomposition, e„ — t- 0, and the stated behaviour of the 
Lebesgue measure is valid in this limit. Using these quantities, labelled by the cells of a cell 
decomposition, we can smear any scalar density of weight {1 — w) and construct a function on 
the configuration/phase space. 

For a Grassmann variable ^ (the spinorial index is suppressed) of weight 1/2, we define the 
Grassmann valued functions 

Under a diffeomorphism: i?„ — )• B'^, its Lebesgue measure — ^ = l^kn) d^^' = \ ^\d^x 
while ^'{x') = J\^\^{x). This implies that 



/3/2 / 3/2|a^.'M/2 



3 Xn{x)^{x) _ 
ax - tJn- 

Thus, the 9 variables are indeed invariant under diffeomorphism^*^. It is immediate that 9n = 

{en)*. 

In the limit of finer cell decomposition, en — )• 0, we get the following expressions: 
lim Xe„ivn,x) = 5y^^x (Kroneker delta), 
lim ^^"^''"''^^ := S'^ivn^x) (5-distribution). 

From these it follows that 

{emA}=j d'x [ d^^^mx)My)} 

f 3 z^^^"^) 

— i^m,n / d X 2 — i^m,ni ■ Xm{,x)Xn{,X^ — ^mnX.m{x)- 

With these, we obtain the Poisson brackets among the Q variables (functions on the phase space). 
The identity in the third line above holds because the Xn^x) = unless x is in the interior of 
the cell Bn- 

At finite but small e^, we can write, Bn ~ ul'^'^ivnj. Notice that the 0„ variables are 
effectively associated with the points Vn- This is particularly useful in expressing the constraint 
expressions in terms of the 0's as illustrated below. 

Recall from equation (3.3) that the smeared Hamiltonian constraint has fermion bilinears 
with and without derivative as well as terms quartic in fermions. Using a cell decomposition, 
we can express the integral as a sum by restricting the integrals over the cells. In the limit of 
small e^, these integrals can be approximated by using the mean value theorem. The factors 



^°In the units where ft = 1 = c, these variables are also dimensionless. This is because the action is dimensionless 
which implies ^' ~ L~^^^ and e ~ L. This is similar to the holonomy variables constructed from connections. For 
the similar variables corresponding to scalar fields, this is not so. 



Matter in Loop Quantum Gravity 



17 



of e^'s available can be distributed with the fermions to go from to 9n variables. The 

correct density weight of +1 ensures that no factors of en's remain unabsorbed. In equations 



n 

n n 
n 

This takes care of the fermionic bilinears. The quartic terms have a factor of y/q{x) in the 
denominator. Following [6, 14, 15], we first express it as a Poisson bracket with only positive 
powers of the volume of the cell containing the point x. From the definition of the determinant 
and our conventions given in the footnote 4, it follows that 



sgn(e)_l,,, 



Let Vn ■= d^Xy/q{x). Then 

1 f ^T/«V"^/^ if T P R 

{AU-)M/'} = \v-nK{x)M = 6 ^'^^^ ' ^ 

1^0, otherwise , 



For small e„, and for non-zero answer (x G Vn ~ ^^n^/^i^) ^'^^ therefore V^/g(x) e^. 

These factors of are neatly combined with the 4 ^ variables to go over to the 6 variables. 
The from the integration measure combines with the three A\ to produce the combinations 
Tt{£^'^^hJ^{hi, . . . }hj^{hj, . . . }h]^^{hK, ■ ■ • }) with the £ ensuring that the three edges are 
non-coplanar and the Tr ensuring the SU(2) invariance. 

Thus all integral expression in the constraints are expressible in terms of the 9 variables. 
Since the transition from the ^ variables to the variables is a linear operation, the Grassmann 
properties of the 9 variables are the same as those of the ^ variables. We have thus the classical 
phase space for Grassmann valued fields. 

To construct a quantum theory, we first construct the space of complex valued (and Grass- 
mann valued) functions and define an inner product on it. The elementary 9 variables, are 
to represented as multiplicative operators while the their complex conjugates ("momenta") are 
to be represented by derivative operators. The procedure is standard [13, 22, 58, 59, 60, 61]. 
The Grassmann nature restricts the functions to be polynomials eg for a single pair of complex 
Grassmann variables, 9, 9, the most general function is: f{9, 9) = a + b9 + c9 + d99 where a, 6, c, 
d are complex numbers. A holomorphic function, has no 9 dependence. For k number of complex 
pairs, the vector space of holomorphic functions is 2'^-dimensional. On this space one has the 
usual Berezin measure, J d9d9f{9, 9) = d, however as pointed out by Thiemann, this is not pos- 
itive definite and thus unsuitable for constructing a Hilbert space. The modification proposed 
by Thiemann is: fi{9,9) := e^^d9d9. The inner product defined by {f,g) := J fi{9,9)f*{9)g{9), 
turns the function space into a Hilbert space. 
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On this, define the operators 

[em:=ef{e), [§f]{e) :=h^. 

These operators satisfy: [9,0]+ = h as desired by the general quantization rule, {A, B} — )■ 
{—i/h)[A,B], valid for the generalized Poisson bracket. It follows that 6 = {6)^ also holds 
reflecting the classical relation 6 = 9* and iJ,{9, 9) is the unique normalised measure selected by 
the above adjointness property [58, 59, 60, 61]. Generalization to finitely many Grassmann vari- 
ables is immediate. For fermionic fields (infinitely many variables) one proceeds via projective 
limit. 

For d-Weyl fermion fields, we will have 2d 9 variables, 9^,9"^, ... , 9'^'^ at each point of S3. The 
holomorphic functions are polynomials of maximal degree 2^*^ in the 9 variables at each point. 
Denote a monomial of degree A: at a point v by 

■■=9''{v)9''{v)---9''^(^)iv), h,...,ik G [l,2d], k{v) G [0,22-^]. 

For every finite set of points^^, {v} = {vi,V2, . . ■ , Vn] and a corresponding vector of labels, 
h ■= {^}fc(D2)' • ■ ■ ' {OfcK)}, define the elementary functions of 9 's, 

j=l,n 

Finite linear combinations of these elementary functions are called (fermionic) cylindrical func- 
tions which are orthonormal with respect to the inner product defined point-wise. Closure of the 
set of these cylindrical functions defines the Hilbert space of the loop quantization of fermions. 
Quantization of more general observables proceeds by first expressing these observables in terms 
of the 9 variables and promoting them to the corresponding operators. 



3.5 Implications for loop quantum cosmology 

Implications of fermionic matter has also been discussed in the context of homogeneous models 
by Bojowald and Das in [14, 15]. The symmetric fermionic fields are only restricted to be 
constants on the spatial manifold by homogeneity. Isotropy however requires these constants to 
be zero. One way to see this is to note that from the presumed non-zero ^''s we can construct 
spatial vectors as V^^^'^'^ and by isotropy there cannot be any non-zero constant vector. For 
homogeneous, anisotropic models however non-zero fermionic constants are allowed^^. 

This has a major implication for diagonal models. Recall that the homogeneous connection 
and triad are of the form: ^^(x) = <I>}ci;^(x), Ef{x) = yfgQPfXj[x) where, are the Maurer- 
Cartan forms and their dual vector fields. The sub-class of diagonal models is defined by the 
restriction: <!)} := c/A}, P/ := p^Af (no sum over I), where AfAj = (5j, K\K^jk\-£ijk = Sijk- 
With this choice, the homogeneous, gravitational part of the Gauss Law constraint vanishes 
identically. With the possibility of non-zero, homogeneous fermionic degrees of freedom, the 
homogeneous Qp 7^ and therefore the gravitational part cannot be zero. Consequently, the 
connection and the triad variables cannot be diagonalised simultaneously. 

A way out, suggested in [14, 15] is to use different A's for connection and triad. The corre- 
sponding dynamical system is still complicated and a specific choice is made to get a simplified 

^^When matter and gravity are considered together, it is natural (sufficient) to take the set of points to be the 
vertices of graph labeUing a spin-network function in the gravity sector (and matter gauge sector if present). 

^^Only fermionic matter is sensitive to the action of the gravitational gauge group (subgroup of local Lorentz 
transformations). All other matter fields are scalars with respect to this action and thus do not contribute to the 
Gauss constraint. 
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model, wherein the constant fermionic vector is taken to be along, 1-axis. This permits a iden- 
tification of canonical coordinates, 4 in number. To explore parity invariance more explicitly, 
the model is further specialized to Bianchi I LRS class. The main conclusion drawn is that in 
this specific case, parity invariance at the level of underlying difference equation is violated only 
if there are parity violating interactions in the matter sector. However, the difference equation 
remains deterministic even when parity is violated [14, 15]. 



4 Gauge fields 

The construction of the kinematical Hilbert space for gauge fields proceeds exactly as in the case 
of the gravitational, SU(2) gauge fields. For identical logic, in a background independent context, 
the gauge fields are quantized through the use of the holonomies and fluxes of the electric fields. 
The new feature is the structure of the corresponding contributions to the Hamiltonian constraint 
and of course the new Gauss constraints corresponding to the additional gauge invariances. As 
remarked in the footnote 12, there is no contribution to the rotation constraint. The contribution 
to the diffeomorphism constraint is of the same form as in the gravitational sector. The couplings 
of these matter fields to the gravitational ones are through the triad variables for minimally 
coupled matter and there are also non-gravitational couplings among the matter fields. The 
underlying background independence which forces the use of holonomy and fluxes impacts the 
quantization of these interaction terms as well. Loop quantization of gauge fields has been 
discussed by Thiemann [58, 59, 60, 61] and Ashtekar-Lewandowski [6] (see also [1]). The salient 
points are noted below. 



4.1 Maxwell field theory 

The gauge group in the Maxwell case case is Abelian which simplifies some of the technical 
details. In particular, since the holonomies are simple numbers (not matrices), the Poisson 
bracket of holonomy with 3-dimensionally smeared electric field is proportional to the holonomy 
and the 3-smeared electric field maps cylindrical functions to cylindrical functions [6]. We will 
however use 2-dimensional smearing and the fluxes Poisson commute in both cases. 

The classical analysis beginning with the Maxwell action followed by canonical formulation 
using the time gauge parametrization is summarised below 



H(N,N'',A)= [ — ^(P"P'' + B"B^) + — F„bP'' + — 5aP", 
y Svr ^ ^ Svr Svr ' 



8vr^ 

:= ^f^^Tfec, E'' := P° = N F^" 

with Aa, P'' being the canonical coordinates. Note that P", B" are both of density weight 1. As 
usual, before quantization, we need to express the constraints (and other observables of interest), 
in terms of the holonomies and fluxes. 

We have the common factor of qab/\/Q- As before, introduce a cell decomposition with 
cells Bn, containing points Vn and with Lebesgue measure and define Vn ■= d^x^J~q{x). 
Noting that 

{K{^)yn] = lVl-'^sgn{e)V:, I G (0, 1), V„ « ei^{vn) for small e„, 
we write 

^(^) ^ j^elS,,{Al{x),^}{Ai{x),y^}. 
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Next, observe that for 2-smeared fluxes of P", we have #5 := fg P"'£abcdS^'^ P"(f)na(5)e| 
where S" is a (small) surface, u is a point in it and e| (small) is its Lebesgue measure. The 'normal' 
to S is given by, na{S) := £abc^-^^£"^ where a = 1,2 denote the local coordinates on S 
while x"(^) denote the embedding of the 2-surface in S3. For a small triangular surface, we can 
also replace na(5')e| ^£abc^Sj6Sj where 6Sj, SSj denote coordinate lengths of a pair of edges 
of the triangle. Likewise, for an infinitesimal curve e, he ■= Pexp(J A) = exp (J A) implies 
h-^he, . . . } « {/^ A, . . . } « {6te'^Ain, ...}. 

To regularise the Hamiltonian, we introduce a cell decomposition and approximate the in- 
tegral by a sum-over-cells. This provides a and the qab/\/Q expression provides another 
factor of for each term. The 6 factors of e can be distributed as 2 each for the two P's 
and 1 each for the two gravitational connections in the Poisson brackets. To express the vari- 
ables in terms of fluxes and holonomies, we need to choose two surfaces and two curves paying 
attention to the contraction of a, 6 indices which is part of specification of gravitational cou- 
pling. 

We recall from the construction done in the gravitational sector. We choose a triangulation 
of S3 by elementary tetrahedra A„ [58, 59, 60, 61]. In anticipation of the quantization step, we 
choose the edges of the tetrahedra to be analytic. One vertex of each tetrahedron is distinguished 
and at the most finitely many tetrahedra meet at such a vertex. Edges meeting at a vertex are 
taken to be out-going. The three edges of each tetrahedron, provide three linearly independent 
tangent vectors at its vertex e°'j, I = 1,2,3, such that the a^^ coordinate interval of the 
I^^ edge = dtje'}, no sum over / and / is not to be confused with the Lorentz index which 
is no longer relevant now. 6tj refers parametrization of the edge. Each tetrahedron, also 
provides three non-coplanar surfaces Sjj, bounded by the edges e/, ej and the 'opposite' edge. 
We have thus a natural choice of surfaces and edges. For these surfaces and edges, we have, 
$/j « lP''£abc{e'iStj){e'jStj) and h]^\hK,VV^} ~ {StKe^Aln, Clearly 

^£'^''<^ijh],'{hK,v^} « 3^^"'^P'^^afec(5^/e5)(,5^JeS)((5^x4){^>^,\/K}• 

Now 

£^-^^e\e''jei = £^'"^ det(e?), det(e?)(5ti5t25t3 =: e^ £abc£^'"^ = 26^ 
Using these, it follows 

I X3 ^7=^^^^' " ^ M^'^n^^ v^}) v^}) 

/ 4 \ glJK gl'J'K' 

X Tr [/i^i { Hk , ^A^} h]^] { hK' , Tvl} . 

In the last equation, we have included TiTiTj = All the factors of e's have been 

neatly absorbed, again thanks to density weight 1 of the Hamiltonian. The sum over n is 
a sum over vertices of the triangulation. The sum over A„ includes those tetrahedra whose 
distinguished vertex is Vn and this is a finite sum. Note that some vertices of the triangulation 
may not be distinguished vertices of any tetrahedron. Such vertices do not contribute to the sum. 
Therefore, the sums over Vn and A„ together, denote sum over the cells of the decomposition. 
Notice that surfaces over which the electric flux is taken is correlated with the edge along which 
the gravitational holonomy is to be taken thanks to the £^-^^ factors. 
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In the passage to quantum theory, the triangulation is adapted to the underlying graph of 
a cyHndrical function^^ by ahgning the edges of the elementary tetrahedra to be segments of 
corresponding edge of the graph and the distinguished vertex to be the vertex of the graph 
[58, 59, 60, 61]. Upon quantization, only those tetrahedra will contribute which have an edge 
overlapping with the edge of the holonomy. In each contributing cell (tetrahedron) , the action of 
the electric flux operator, #/j will be non-zero only on the Maxwell holonomy along the edge ex 
which is transversal to the surface Sjj and will give hqx- This also selects the edge K in the 
gravitational holonomy, hx- Here qx, an integer, is the charge (in units of electric charge) in the 
Maxwell holonomy. On any cylindrical function, the action of the regulated quantum operator, 
results in a finite sum of terms since only vertices of the underlying graph contribute to the 
sum, regardless of refinement of the triangulation. 

We have deviated somewhat from the procedure given in [6, 58, 59, 60, 61]. We have not done 
the point splitting and explicitly made passage to the flux operators absorbing the factors of e's 
(in the limit of small e) . The correlation between the edge in the gravitational holonomy and the 
surface in the Maxwell flux is preserved through the use of triplet of edges and corresponding 
unique transversal surfaces provided by an elementary tetrahedron. 

For the magnetic term, we can proceed identically by defining a flux *5 := B°'£abcdS'"^ 
which will lead to an equation same as the one above with #5 replaced by ^5. Unlike the electric 
flux however, this magnetic flux is not an elementary variable. For 



where es is a small close curve bounding the surface S. When the triangulation is adapted to the 
underlying graph of a cylindrical function, the curve es may be taken to be one of the loops ajj 
of Thiemann [58, 59, 60, 61], which bounds the face Sij. Because of the presence of the Maxwell 
holonomies (with unit charge), the underlying graph of a cylindrical state is changed by adding 
these loops. But there are no gravitational holonomies along these loops and hence the spin 
labels of these extra loops are zero. 

We note in passing that if we have a spin and charge network state, i.e. a graph with both 
the spin and the charge labels for its edges some of which could be zero, then the action of the 
electric part of the Hamiltonian will be non-zero only on those edges which have both the spin 
and the charge labels to be non-zero - the volume operator will kill edges with no spin while 
the electric flux operator will kill the edge with no charge. This feature enables the operator 
to be well defined on the kinematical Hilbert space even in the limit of infinite refinement. By 
contrast, the magnetic part will not kill an edge with zero charge but will add loops with unit 
charge. This prevents the operator to be well defined on the kinematical Hilbert space in the 
limit of infinite refinement. It could however be defined on a subspace of the algebraic dual of 
the space of the cylindrical functions [6]. 

4.2 Yang-Mills field theory 

Apart from the non-Abelian nature of the gauge field, in form, the expressions proceed in the 
same manner as for the Maxwell case. The details are given in [58, 59, 60, 61]. 

^^The cylindrical functions in this case are analogous to the gravitational spin network functions except for the 
replacement of the spin label by a integer charge label with charge conservation at each vertex and are called 
charge network functions [62]. These have been called flux networks in [6]. 
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5 Propagation of matter waves on quantum geometry 

While there are good arguments to justify search for a quantum gravity, its observable signatures 
have been quite difficult to come by. Various different approaches to quantum gravity have 
contained a suggestion that space-time of quantum gravity is different from the continuum 
classical geometry be it discreteness and/or non-commutativity etc. Coupled with a heuristic 
idea that existence of a fundamental length scale may indicate a conflict with Lorentz invariance, 
a possible violation of Lorentz invariance has been considered as a possible signature of quantum 
gravity. The violation is thought to be manifested as a deviation from the Lorentz invariant 
dispersion relation for particle propagation on a 'quantum space-time'. Its cumulative effect 
over cosmological distances are estimated to be detectable. From LQG Hamiltonian constraint, 
such deviation have been extracted [1, 2, 3, 26] using the following set of assumptions^^. The 
Lorentz violations in the Polymer quantized scalar in Minkowski background has already been 
discussed in Section 2.5. 

One begins by assuming a quantum state to be of a product form with a suitable state 
for geometry as well as for matter. For example, the early works, e.g. [26], the geometry 
state was a weave state [10, 17, 28, 36, 37, 38, 65] which provides a length scale Lweave S> 
such that for larger scales the 3-geometry may be approximated by a continuum while the 
discrete structure manifests near the Planck scale. The scales relevant for observations of interest 
are assumed to be much larger than L„eave- The matter state is chosen to be a coherent 
state so that the expectation values of the matter operators can be approximated by classical 
fields (fields obeying classical equations to leading order). These fields are assumed to be 
varying slowly over a length scale of A ^ L„eavc- Taking expectation value of the Hamiltonian 
constraint, the geometry operators are replaced by the metric used in the construction of the 
weave plus corrections of the order of ip/L^eave- The expectation values of the matter part of 
the Hamiltonian constraint is taken as defining an effective Hamiltonian. The space-time metric 
is taken to be a static metric determined by the 3- metric and a lapse A^. The matter equations of 
motion following from the effective Hamiltonian contain the modifications implied by quantum 
gravity. In computing the effective Hamiltonian, one expands the classical fields at the vertices in 
a region of size Lweave around a central point of the region and averages over these vertices. The 
symmetries of the 3-geometry (eg rotational invariance of the flat geometry) is used to restrict the 
form of the averages. In turns out that there are corrections in the quadratic order in the Maxwell 
fields which imply polarization dependent speed of propagation suggesting birefringence. The 
computations are varied by different choice of states for expectation values eg more general 
semi-classical states rather than weaves or even using only qualitative properties of 'would be 
semi-classical states' [1]. These computations also give additional corrections. Subsequently, the 
method is also applied to propagation of neutrinos [2, 3]. These early computations show how 
Lorentz invariance violations - as manifested by dispersion relations - could potentially arise in 
LQG. But do they imply that LQG violates Lorentz invariance? 

This is not easy to answer. The above computations are essentially at a kinematical level. 
This however need not be regarded as a strong drawback, because the very characterization of 
Lorentz violation pre-supposes a background space-time in which one can find locally inertial 
observers. Introducing a background is akin to a gauge fixing and heuristically at least the use of 
kinematical states could be justified. There are still not satisfactory semi-classical states of the 
geometry and matter system which will describe 'test fields' on a given space-time background. 
The proper framework for this would be quantum field theory on curved space-times. To obtain 
it from the full quantum theory of gravity and matter, is a challenging task. First steps in this 



The subject of Lorentz invariance violation in all its aspects is rather large. We have included it here only to 
the extent of indicating how a computation of potentially observable signatures involving matter could be done. 
For review, we refer the reader to [46] . 
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direction are discussed by Sahlmann and Thiemann [55, 56]. A cautionary note on a different 
aspect, has also been sounded in [16], pointing out that staying within Hamiltonian framework 
for computation including approximation leading to higher order derivatives, could be misleading 
because the corrections could be due to additional degrees of freedom with a Lorentz invariant 
dynamics. There is yet another aspect. Dispersion relation, relativistic or otherwise, need 
a background space-time. The presumed space-times in the first approaches have been static 
space-times with spatial geometry being described by a semi-classical state. Thus the 'time' part 
of the geometry does not have anything to do with the quantised spatial geometry. One has 
not used a quantum space-time, only quantum space. An extension to quantum space-time has 
been considered in the context of FLRW cosmology [5] . Such extensions are probably needed to 
discuss matter wave propagations on quantum space-time and its implications. 

6 Concluding remarks 

We began by noting that non-Fock quantization of matter fields is necessary once 'space-time' 
is dynamical. However, the illusion of a fixed background space-time with Fock quantized 
fields on it, is also extremely persuasive. In these circumstances, gravitational interactions of 
matter fields are negligible compared to intra-matter interactions. Consequently, the observed 
semi-classicality of geometry itself would conceivably be sensitive to matter content as well as 
their interactions. Matter interactions must play a role in getting our 'late universe' with its 
essentially classical space-time emerge from the quantum gravity regime of the early universe. 
Study of (loop) quantized matter has the potential to provide selection criteria from the stability 
of semi-classical states of geometry. We are still quite far away from this goal. 

There are other 'applications' of matter fields. They can serve as probes for understanding 
properties of background independent quantization. For instance, in a fixed background space- 
time such as the Minkowski space-time, the short distance behaviour of quantum fields play a role 
in admitting the chiral (and the trace) anomalies. One of these, namely the axial anomaly, has 
played a historic role in understanding the vr" — )• 27 process. Cancellation of gauge anomalies 
has also provided constraints on model building. In a background independent quantization, we 
loose the handle of short distance behaviour. How then are the anomalies to be understood? This 
is an old issue more recently mentioned by Nicolai et al. [52], still awaiting an understanding. 

The focus of analysis of matter in quantum gravity has been to study its propagation in the 
regime of semi-classical gravity. These have indicated that one could expect dispersion rela- 
tions not compliant with Lorentz invariance. This has been seen in the propagation of massless 
fermions (neutrinos), light as well as massless scalar fields. For fermions and light, quantum 
effects of spatial geometry are incorporated while space-time is constructed classically. The 
case of the scalar field discussed in Section 2.5, is qualitatively different from these. Here the 
relativistic (massless) scalar field is classically expressed as a collection of classical 'harmonic 
oscillators' with frequencies Wk ~ |k| and there is no Lorentz violation at this stage. For the 
usual Fock quantization, the energy spectrum, in units of Hu^ is linearly spaced and the Lorentz 
invariant dispersion results from the propagator. With polymer quantization, the spectrum is 
not linearly spaced and leads to Lorentz violating dispersion relation. While one can construct 
Lorentz generators from the VTk, (/>k basic variables, one cannot do so with the C/Ak — t^^k '^^^'i" 
ables for A 7^ 0. Clearly, it is the polymer quantization which requires the use of C/^k with 
non-zero A which is the source of this Lorentz violation at the level of 'free field'. Importing the 
Lorentz violating dispersion relation into the usual computational scheme of perturbative quan- 
tum field theory based on Fock quantization with its UV divergences (especially the 'quadratic 
divergences'), will lead to the 'fine tuning' problem discussed by Collins et al. [18] (see also [40]). 
To check the viability of a polymer quantized field theory on a Minkowski background, one will 
need to develop a computational scheme for a polymer quantized, interacting theory and then 
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test if at 'low energies' Lorentz violating effects are suppressed or not. Currently this lies only 
in the realm of possibilities. 
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